Vortex-enhanced alternating order around impurities in antiferromagnets 
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It has been recently pointed out that the presence of a nonmagnetic impurity in a Heisenberg 
antiferromagnet generates an alternating order of the surrounding spins which is independent of 
temperature in a wide range. Quantum Monte Carlo simulations in the two-dimensional 5 = 1/2 
case confirmed this picture, but showed a counterintuitive enhancement of the alternating order 
around the Kosterlitz-Thouless transition. We propose here an explanation in terms of the effect of 
vortex excitations. 



In a study of the effects of impurities in Heisenberg 
antiferromagnets, Eggert, Syljuasen, Anfuso, and An- 
dres [l[ used an effective model in order to derive a 
universal expression for the function (n z (r)) describing 
the local alternating magnetization induced around a va- 
cancy at site r — by a uniform magnetic field B applied 
along the z-axis. Basically, (n z (r)) behaves as e~ r l x 
and the main point of interest lies in the fact that the 
characteristic length A(T) is actually independent of 
the temperature T, as they also checked in a wide tem- 
perature range by quantum Monte Carlo simulations of 
the two-dimensional S=\ Heisenberg antiferromagnet. 
In addition, they point out an 'exotic effect' observed 
in their simulations: when T steps over the estimated 
Kosterlitz-Thouless temperature Tkt the temperature 
behavior of the induced alternating order is nonmono- 
tonic, i.e., A(T) first increases and displays a little bump 
before being progressively reduced as non-linear thermal 
fluctuations rapidly destroy long-range correlations. The 
last sentence of Ref. [l| claims that this feature is 'coun- 
terintuitive and calls for further investigation'. Here we 
aim indeed at giving an explanation in terms of the ef- 
fects of vortex excitations. 

In what follows all the symbols, if not explicitly de- 
fined, will have the same meaning as in Ref. [l[ . 

The ingredients of the effective model of Ref. [l[ are 
the spin stiffness p s and what we here call the effective 
elastic constant G = B 2 /4zJ = G u . The latter is obtained 
by studying the stability of the uniform spin-flop config- 
uration, and it comes out that A = y/ p s /G. 

However, if we allow for a 'vortex' configuration the in- 
plane components are nonuniform, their orientation be- 
ing described by a local azimuthal angle <pi. By following 
the same path of Ref. |lj , after setting 

( s iT s i)— ± s cos(a±<5) (cos</3i,sin(/?j) 

and sf = s sm(a±S) for the two sublattices (±), the 
single- vortex correction to the effective energy 

5E e s = /O (cos 2a + cos 28)C , 

is obtained Q; in the last equation we have set 



JC = z.Js 2 /2 and 



C = ^TnY. [1 - cosfo - ¥>i)] 
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where the second line follows from using the continuum 
limit for a single vortex centered in ro, so that |Vp(r)| = 
1/| r — 7*o | ; ttR 2 represents the lattice surface available 
for the free-vortex. Therefore, the full effective energy 
becomes 



-JC(1 - 2C) - 2K,C sin 2 a + 2/0(1- C) sm 2 8 
—Bs cos a sin (5 : 



minimization with respect to the canting angle S gives 
4/0(1 — C) sin S — Bs cos a, and the resulting effective en- 
ergy acquires the same form obtained in Ref. [l| but for 
the effective elastic constant that now reads: 
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2/0/3 = G v 



with a relative deviation from the uniform case given by 
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For T < Tkt vortices are tightly bound in pairs and 
i?~l, so that G v — G u . At Tkt the first vortex- 
antivortex pair unbinds and R ~ ^/N. A rapid devia- 
tion of G from the uniform value G u appears and the 
effect reported in Ref. [l[ can be reasonably observed. 
Indeed, the parameter values used in the inset of Fig. 2 of 
Ref. [l|, 7V = 128 2 and S/J = 0.2, give £~3.8-lCT 4 and 
2z Js/ B ~ 20, resulting in 8G/G~ —0.15 and eventually 
8\/\ ~ + 0.075, a variation that is compatible with what 
was observed in [l[ . The initial increase of the character- 
istic length A at Tkt can thus be reasonably traced back 
to free-vortices softening the elastic force that induces 
the spin-flop configuration. 
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A more quantitative description of this effect is ap- 
parently complicated and beyond the scope of this note. 
Nevertheless, one can figure out that the exponentially 
fast decrease of the in-plane correlation length £ (T) above 



Tkt would very soon impose an upper bound to A in- 
ducing its progressive decrease as temperature is further 
raised. 
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in a dilute gas approximation the respective contributions 
simply add, and C — ► nC, n being the total number of 
vortex-like excitations. 



